David  Taylor  Research  Center 

Bethesda,  Maryland  20084-5000 


AD-A241  208 


DTRC  /  SHD-1 350-01  March  1 99 1 

Ship  Hydromechanics  Department 


A  MATHEMATICAL  MODEL  FOR  SURFACE  SHIP  MANEUVERING 

by 

William  R.  McCreight 


91-11291 

Illlllltlllllilllllll 


Approved  for  Public  Release 
Distribution  is  Unlimited 


MAJOR  DTRC  TECHNICAL  COMPONENTS 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  PLANS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1 .  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


.  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 

-a  REPORT  SECURITY  CLASSIFICATION 

lb.  RESTRICTIVE  MARKINGS 

Unclassified 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  Public  Release 

2b  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

Distribution  is  Unlimited 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

DTRC  /  SHD-1350-01 

£a  NAME  OF  PERFORMING  ORGAN, ZATION 

David  Taylor  Research  Center 


6b.  OFFICE  SYMBOL 
(If  appicable) 

1562 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Naval  Sea  Systems  Command 


6c.  ADDRESS  (City,  State,  and  Zip  Code) 

7b.  ADDRESS  (City,  Stale,  and  Zip  Code) 

Bethesda,  MD  20084-5000 

Washington,  D.C.  20362 

8a.  NAME  OF  FUNDING  /  SPONSORING 

8b.  OFFICE  SYMBOL 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

ORGANIZATION 

(tt  applicable) 

Naval  Sea  Systems  Comm.md 

SEA  55W3 

8c.  ADDRESS  (City,  State,  and  Zip  Code) 

Washington,  D.C.  20362 


1 1  TITLE  (Include  Security  Classification ) 

A  Mathematical  Model  for  Surface  Ship  Maneuvering 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 
243 13N 


PROJECT 

TASK 

NO. 

i _ _ _ 

NO. 

WORK  UNIT 
ACCESSION  NO 
DN5G8317 


12.  PERSONAL  AUTHOR(S) 

William  R.  McCreight 


13a  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 


V  I 


14.  DATE  OF  REPORT  (Year,  Month.  Day)  I  1 5.  PAGE  COUNT 

1991  March  I  24 


COSATI  CODES 


GROUP 


SUB-GROUP 


1 8.  SUBJECT  TERMS  (Continue  on  reverse  II  necessary  and  Identity  by  block  number) 

Maneuvering 


1 9  ABSTRACT  (Continue  on  reverse  if  necessary  and  Identll )  by  block  number) 

The  mathematical  model  for  a  time  domain  simulation  of  a  surface  ship  maneuvering  in  calm  water  is 
presented.  The  six-degree-of-freedom  mathematical  model  is  applicable  to  conventional  monohulls  or 
SWATHs.  Included  are  calm  water  hydrodynamic  forces,  hydrostatic  forces,  unsteady  wind  forces, 
slowly-varying  wave  drift  forces,  and  forces  due  to  a  towed  body.  The  model  depends  upon  data  derived 
from  model  experiments.  This  model  represents  an  improvement  over  previous  models  in  several  respects. 
Most  notable  are  the  changes  in  the  modeling  of  resistance,  calm  water  surge  hydrodynamic  forces,  and 
propeller-rudder  interaction  which  together  lead  to  improved  speed  loss  prediction. 


20  DISTRIBUTION  /  AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIED  O/UNUMITED  0  SAME  AS  HPT 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

William  R.  McCreight  _ _ 


DD  FORM  1473,  *4  MAR 


D  OTIC  USERS 


21 .  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22c.  OFFICE  SYMBOL 

1561 


13  APR  adWon  may  ba  un4  untH  axteualad 
Ai  otter  adrttona  ara  obaotota 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


i 


UNCLASSIFIED _ 

SECURiTV  CLASSIFICATION  OF  THIS  PAGE 


SECURITY  CLASSIFICATION  OF  THIS  PAGF 


11 


CONTENTS 


page 


ABSTRACT . 1 

ADMINISTRATIVE  INFORMATION . 1 

INTRODUCTION . 1 

EQUATIONS  OF  MOTION . 4 

COORDINATE  SYSTEMS . 4 

DYNAMICS . 4 

EXTERNAL  FORCES  AND  MOMENTS . 6 

HYDROSTATIC  FORCES . 6 

CALM  WATER  HYDRODYNAMIC  FORCES . 7 

Resistance  and  Propulsion  Model . 9 

Machinery  Model . 11 

WIND  FORCES . 12 

Unsteady  Wind  Model . 12 

Approximation  for  Wind  Force  Coefficients . 13 

WAVE  FORCES . 14 

Slowly-Varying  Wave  Drift  Force  Model . 16 

TOW  FORCES . 17 

SOLUTION  OF  THE  EQUATIONS  OF  MOTION . 18 

CONCLUSION . 19 

APPENDIX  -  ADDITIONAL  TERMS  IN  WAVMAN44  OR  TOWMAN87 . 21 

REFERENCES . 23 


in 


ABSTRACT 


The  mathematical  model  for  a  time  domain  simulation  of  a  surface 
ship  maneuvering  in  calm  water  is  presented.  The  six-degree-of- 
freedom  mathematical  model  is  applicable  to  conventional 
monohulls  or  SWATHs.  Included  are  calm  water  hydrodynamic 
forces,  hydrostatic  forces,  unsteady  wind  forces,  slowly-varying  wave 
drift  forces,  and  forces  due  to  a  towed  body.  The  model  depends  upon 
data  derived  from  model  experiments.  This  model  represents  an 
improvement  over  previous  models  in  several  respects.  Most  notable 
are  the  changes  in  the  modeling  of  resistance,  calm  water  surge 
hydrodynamic  forces,  and  propeller-rudder  interaction  which 
together  lead  to  improved  speed  loss  prediction. 

ADMINISTRATIVE  INFORMATION 

this  report  was  funded  by  Program  Element  24313N,  D  N  Number  DN  508317  and  Work 
Unii  Number  1-1235-985 


INTRODUCTION 

This  report  documents  the  current  state  of  a  time  domain  simulation  for  the 
maneuvering  of  surface  ships,  including  the  effect  of  unsteady  wind  forces  (wind  gusts),  slowly- 
varying  wave  drift  forces  in  long-  or  short-crested  seas,  and  forces  due  to  a  towed  body.  This 
model  is  suitable  for  trackkeeping  studies  as  well  as  traditional  maneuverability  studies.  It 
has  evolved  at  the  David  Taylor  Research  Center  (DTRC)  over  a  period  of  more  than  two 
decades  from  the  original  three-degree-of-freedom  model  of  Strom-Tejsen[l],  based  on  the 
mathematical  model  of  Abkowitzl2],  Many  people  have  worked  on  it,  or  influenced  its 
development. 

Motter1  added  a  fourth  degree  of  freedom  (roll)  and  the  propeller  and  rudder  model  of 
Smit  and  Chislett  [3].  McCreight  and  Motter1  carried  out  a  major  rewrite  of  this  program  as  the 
first  stage  of  an  effort  to  develop  a  full  maneuvering  in  waves  model.  In  addition  to  structuring 
the  program,  a  steady  wind  model  was  added,  coefficients  were  added  to  be  valid  over  an 
extended  speed  range  instead  of  merely  accounting  for  "small"  deviations  from  the  initial 
conditions,  improved  roll  damping.  Provision  was  made  for  defining  an  arbitrary  sequence  of 
rudder  and  speed  commands.  A  nonlinear  initialization  scheme  was  developed  to  handle  the 
large  excursions  which  can  be  caused  by  strong  winds.  Finally,  the  previously  used  Euler 
integration  was  replaced  by  an  adaptive  Runge-Kutta  which  allowed  the  use  of  large  step  sizes 


1  Published  in  reports  with  limited  distribution. 
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while  maintaining  accuracy.  This  was  particularly  important  due  to  the  sensitivity  of  the 
model  in  roll.  In  addition,  a  detailed  model  of  steam  turbine,  gas  turbine,  and  diesel  engines 
was  incorporated  by  Propulsion  Dynamics,  Inc.  (PDI)[4).  PDI  also  produced  a  user's  manual  for 
this  program.  This  maneuvering  program  was  known  as  WAVMAN44  for  historical  reasons, 
even  though  it  has  no  wave  effects.  WAVMAN44  has  been  used  in  many  maneuvering  studies, 
and  has  been  the  base  for  two  further  developments. 

The  first  and  most  extensive  of  the  developments  is  the  originally  planned 
maneuvering  in  waves  simulation,  which  added  two  additional  degrees  of  freedom  (heave  and 
pitch)  to  produce  a  full  six-degrees-of-freedom  and  added  first  order  wave  excitations  and 
reaction  forces,  including  memory  effects  (which  result  in  the  frequency-dependent  added  mass 
and  damping  effects  in  the  frequency  domain).  Regular  wave,  and  long  and  short-crested 
random  seas  models  were  included.  The  results  of  this  effort  are  described  by  McCreight[5). 
The  resultant  computer  program  has  been  referred  to  as  WAVMAN. 

The  other  branch  of  development  for  WAVMAN44  was  by  Waters,  Hickok,  Turner, 
Hart,  and  Bochinski  of  DTRC.  They  implemented  a  series  of  efforts  which  are  either  presented 
in  reports  of  limited  distribution  or  are  undocumented.  Their  calm  water  hydrodynamic  model 
was  believed  to  be  more  appropriate  for  SWATH  ships  with  dihedral  fins  forward  of  the 
propellers,  and  included  the  addition  of  fixed  fin  effects  (together  with  heave  and  pitch  for 
examining  the  sinkage  and  trim  changes  due  to  fins),  an  unsteady  wind  model,  a  mean  wave 
drift  force  model,  forces  due  to  a  towed  body,  and  a  simple  autopilot.  Equations  for  calm  water 
forces  for  heave  and  pitch  were  added.  This  program  has  been  referred  to  as  TOWMAN87. 

The  present  version  of  the  program,  MPSS  (Maneuvering  Program  for  Surface  Ships),  is 
applicable  to  both  monohull  and  SWATH  designs.  While  TOWMAN87  served  as  the  starting 
point,  numerous  features  from  WAVMAN  have  been  incorporated  in  MPSS.  MPSS  also  has 
components  and  characteristics  not  included  in  any  of  the  predecessor  programs.  A  list  of  the 
major  differences  and  similarities  between  MPSS  and  the  other  programs  follows. 

(1)  MPSS  is  a  six-degree-of-freedom  model  as  is  WAVMAN.  While  TOWMAN87 
includes  some  heave  and  pitch  calm  water  hydrodynamic  and  hydrostatic  terms,  it  does  not 
include  the  necessary  inertial  terms  and  transformations  between  coordinate  systems  required 
for  a  consistent  six  degrees  of  freedom  model. 

(2)  The  terms  used  in  the  calm  water  hydrodynamic  model  in  MPSS  were  selected  on 
the  basis  of  physical  considerations.  The  terms  in  the  model  were  examined  for  the  symmetry 
or  asymmetry,  resulting  in  changes  in  the  equations.  Fourteen  terms  in  the  TOWMAN87  model 
were  eliminated  and  30  terms  were  added  on  this  basis.  The  model  is  most  similar  to  that  of 
TOWMAN87.  For  flexibility,  the  capability  to  utilize  the  model  implemented  in  WAVMAN 
(and  WAVMAN44)  has  been  retained. 
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(3)  MPSS  incorporates  an  explicit  model  to  account  for  the  effects  of  the  propeller- 
rudder  interaction  on  the  forces  generated  by  the  rudder  similar  to  that  developed  by 
Norrbin[6],  Abkowitz[7],  and  others.  This  is  extremely  important  in  modeling  the  motions  of 
vessels  with  rudder(s)  aft  of  the  propeller(s).  This  configuration,  which  occurs  in  most 
monohulls  as  well  as  SWATHs  with  overhanging  struts,  results  in  the  propeller  wake  flowing 
over  the  rudder.  This  mode!  is  also  applicable  to  ships  with  rudders  forward  of  the  propeller. 

(4)  A  linear  interpolation  method  for  accounting  for  the  ship  resistance  was  adopted  in 

MPSS. 

(5)  Approximations  for  the  wind  forces  from  WAVMAN  are  included  in  MPSS.  The 
unsteady  wind  model  and  the  forces  due  to  a  towed  body  from  TOWMAN87  have  been  retained. 

(6)  A  slowly-varying  wave  drift  force  model  has  been  developed  and  implemented  in 

MPSS. 

(7)  MPSS  has  features  not  present  in  either  WAVMAN  or  TOWMAN87  which  do  not 
affect  the  modeling.  These  include  a  new  data  input  scheme  which  is  more  convenient  and 
verifies  variable  names. 

The  changes  in  the  surge  force  model  in  items  (2),  (3),  and  (4)  above  together  contribute 
tc  a  greatly  improved  speed  loss  model.  Another  effect  of  these  changes  is  that  only  one  set  of 
coefficients  is  required,  which  is  valid  for  all  speeds,  instead  of  separate  sets  at  each  of  the 
speeds  which  is  sometimes  done. 

The  resultant  mathematical  model  which  is  implemented  in  MPSS  is  described  in  this 
report.  As  noted  above,  documentation  of  several  intermediate  stages  of  the  developments  of 
this  model  has  been  in  limited  distribution  reports.  The  distribution  was  limited  due  to  the 
inclusion  of  hydrodynamic  coefficients  and  predicted  maneuvering  performance  for  specific 
designs  in  those  reports.  For  this  reason,  such  detailed  data  has  been  omitted  from  this  report. 
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EQUATIONS  OF  MOTION 

Newton's  laws  of  motion  form  the  basis  of  calculating  the  acceleration  of  the  ship, 
considered  to  be  a  rigid  body  subject  to  external  forces  due  to  the  propulsion  system,  rudder,  fins, 
the  relative  motion  of  the  water  due  to  the  motion  of  the  ship,  and  environmental  forces  due  to 
wind  and  waves.  These  accelerations  are  then  integrated  twice  to  obtain  the  time  evolving 
position  and  orientation  of  the  ship  as  it  maneuvers. 

COORDINATE  SYSTEMS 

In  developing  the  equations  of  motion,  several  related  coordinate  systems  are  used. 
These  are  presented  here  in  the  order  of  the  natural  succession  from  one  to  the  next. 

These  are: 

1)  The  inertial  (xo,  yo,  zo)  system  which  is  right  handed,  with  zq  positive 
downwards,  and  is  fixed  in  space. 

2)  The  upright  (xu,  yU/  zu )  coordinate  system  is  parallel  to  the  inertial  coordinate 
system  with  its  origin  on  the  free  surface,  coinciding  with  amidships  and  moving 
with  the  ship. 

3)  The  yawed  (x’,  y',  z')  coordinate  system  is  obtained  by  rotation  of  the  upright 
coordinate  system  through  the  angle  y  about  the  zu  axis. 

4)  The  pitched  (x",  y",  z")  coordinate  system  is  obtained  by  rotating  the  yawed 
coordinate  system  through  the  angle  0  about  the  y'  axis. 

5)  The  body  or  maneuvering  (x,  y,  z)  coordinate  system  is  obtained  by  rotating  the 
pitched  coordinate  system  through  the  angle  <j>  about  the  x"  axis.  This  is  fixed  and 
is  the  conventional  SNAME  (Society  of  Naval  Architects  and  Marine  Engineers) 
coordinate  system[8]. 

DYNAMICS 

The  equations  of  motion  assuming:  (a)  transverse  symmetry,  (b)  the  principal  axes  of 
inertia  parallel  to  the  chosen  axis  system,  Iy  =  Iz,  and  (c)  the  center  of  gravity  located  at 
(xg,  0,  zq),  are,  in  standard  nomenclature^]: 


X  =  mfii+qw-rv-xG(q2  +  r2)  +  ZG(pr  +  q)]  (1) 

Y  =  m[  v  +  r u-p  w  +  ZG(qr-p ) +  xc(pq  + r )]  (2) 

Z  =  m[w  +  pv-qu-zG(p2  +  q2)  +  XG<pr-q)3  (3) 

K  =  Ixp  +  (Iz-ly)qr-mzG(v  +  ru-pw)-(Ixz  +  mxGZG)(pq+f)  (4) 
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M  =  Iyq  +  (Ix-iz)pr  +  mzc(u  +  qw-rv)-mxG(w  +  pv-qu) 


(5) 

(6) 


+  (Ixz  +  mxczc)  (p2-^) 

N  =  I2r  -  ( Ix  -  Iy  )  pq  +  mxc(  v+ru-pw)  +  (I)c  +  mxczc)(qr-  p) 

In  the  above  equations,  u,  v,  and  w  are  the  surge,  sway  and  heave  velocities,  and  p,  q. 
and  r  are  the  roll,  pitch,  and  yaw  velocities,  measured  in  the  body  coordinate  system.  X,  Y,  Z, 
K,  M,  and  N  are  the  total  external  surge,  sway,  heave,  roll,  pitch,  and  yaw  forces  and  moments, 
respectively,  m  is  the  ship  mass,  Ix,  ly,  and  Iz  are  the  mass  moments  of  inertial  about  the  x,  v, 
and  z  axis,  respectively.  1^  is  the  product  of  inertia  with  respect  to  the  x  and  z  axis. 

The  following  equations  relate  the  time  rate  of  change  of  position  and  orientation  to 
the  velocities  in  maneuvering  coordinates: 

xo  =  u  cosG  cosy  +  v  (  sin<{>  sin9  cosy  -  cos<))  siny  )  +  w  (cosp  sin0  cosy  +  si  no  siny)  (7) 

yo  =  u  cos0  siny  +  v  (  sinO  sin©  siny  +  cos<t>  cosy )  +  w  (cos0  sin0  siny  -  sinO  cosy)  (8) 


z o  -  u  sin0  +  v  sin<(>  cos©  +  w  cos4>  oosfl  (9) 

4>=  p  +  (q  sin<t>  +  "-osO)  tan0  (10) 

0=  q  cose  -  r  jinv  (11) 

y=  (q  sinO  +  r  cosO)  /  cos0  (12) 


When  integrated,  these  equations  ^icld  (xq.  v«,  zo),  the  position  of  the  ship  origin  in  the 
inertial  coordinate  system,  and  the  Euler  angles  0/  6,  and  y  which  are  commonly  known  as  the 
roll,  pitch  and  yaw  angles. 

The  transformation  of  a  vector  from  the  yawed  coordinate  system  to  the  maneuvering 
coordinate  system,  assuming  x  =  z  =  0,  is  given  by 

cos©  0  -sin0 

FmaN  =  sinOsin0  cos4>  sin^cos©  FyAW  (13) 

cos4>sin0  -sinO  cosOcosO 

where  FmaN  and  FyaW  are  the  vectors  in  the  maneuvering  and  yawed  coordinate  systems, 
respectively.  Inis  is  required  io  transform  force  and  moment  vectors  computed  in  the  yawed 
coordinate  system  into  the  maneuvering  coordinate  system. 
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EXTERN  AL  FORCES  AND  MOMENTS 


Combining  the  external  forces  and  moments  into  a  vector: 

F  =  (X,  Y,  Z,  K,  MNt)T 
where  T  indicates  the  transpose. 

We  can  separate  the  forces  into  components  according  to  the  cause  of  each  force  and 
moment: 

F  =  Fhs  +  FCVV  +  FwiND  +  FWAVE  +  FjOW  (14) 

The  subscript  HS  indicates  hydrostatic  force,  and  all  six  components  must  be  considered 
due  to  the  changing  orientation  of  the  ship  axis  system.  The  subscript  CW  indicates  the 
maneuvering  force  affecting  the  hull  in  calm  water.  These  include  the  bare  hull,  propeller, 
rudder,  and  fin,  if  any.  The  next  two  terms  model  external  forces  due  to  environmental 
disturbances.  In  the  present  work,  unsteady  wind  forces  and  slowly-varying  wave  drift  forces 
are  modelled.  The  final  term  models  the  forces  due  to  a  towed  body. 

All  force  and  moment  components  in  the  equation  above  are  each  calculated  in  a 
coordinate  system  which  is  convenient  for  that  component,  then  transformed  into  the 
maneuvering  coordinate  system,  if  necessary,  and  summed  to  obtain  the  total  forces  and  moments 
acting  on  the  ship.  These  models  are  described  in  detail  below. 

HYDROSTATIC  FORCES 

The  hydrostatic  forces  and  moments  which  represent  the  difference  between  the 
gravitational  forces  on  the  ship  and  the  integral  of  the  hydrostatic  pressure  acting  on  the  hull 
are  calculated  using  a  linear  theory  in  the  yawed  coordinate  system.  In  this  system  the  non¬ 
zero  components  of  the  hydrostatic  forces  and  moments  are  taken  as: 

ZHS  =  -20  C33  -  0  C35 


Khs  =  -<>£44 


Mhs  =  -20  C35  -  0  C55 

where  zq  is  the  vertical  displacement  of  the  seakeeping  origin  from  its  rest  position,  and  0  and  0 
are  the  (finite  amplitude)  roll  and  pitch  angles,  which  for  small  angles  reduce  to  the 
corresponding  linear  seakeeping  theory  angles.  The  Qj  terms  arc  the  conventional  linear 
hydrostatic  coefficients  defined  by,  for  example,  Salvesen,  Tuck,  and  Faltinsen[9).  The 
resulting  forces  and  moments  are  then  transforr  led  into  the  body  coordinate  system. 
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CALM  WATER  HYDRODYNAMIC  FORCES 


The  calm  water  terms  include  the  traditional  maneuvering  model  and  have  evolved 
over  many  years.  The  history  and  development  of  these  terms  was  briefly  outlined  in  the 
introduction.  The  calm  water  forces  and  moments  are  represented  as  the  sum  of  the  product  of 
coefficients  and  state  variables.  In  the  computer  program,  all  terms  are  non-dimensionalized 
using  the  SN'AME  system[8].  Following  solution  of  the  system  of  equations  of  motion,  velocities 
are  dimensionalized.  For  simplicity,  the  equations  presented  here  will  be  given  in  the 
dimensional  form.  These  forces  are  calculated  in  the  maneuvering  coordinate  svstem. 

The  final  expression  for  surge,  sway,  and  heave  forces  and  roll,  pitch,  and  vaw  moments 

follow. 


*C\V  Xu^  ^0  ^  Xw  XvvF  F  )  v  +1  ^it  XrrF  F)r-+(  ^vr  Xvvtt  ^  ^  ^  t* 

+  f  XV(t)  V  +  Xj.Q  r  )  <J)  +  (  X^yJ  +  Xyy^p  V^  +  X^Q^,  +  X^QQ  V  T  )  0  ^ 

+  XS8  52  +  XFk5Fk 

+  (  ^cee  +  *xcee£X>  <t>2  )  C2  e2  +  X^  e  0  +  Xrp 
YCW  =  <Yv  +  Yvuu  +YvuuU2)v  +  Yww  +  YpP  +  Yqq  +  (Yr  +  Yruu+Yruuu2)f 

+  (Yv  +  YvFF  +  YvFFF2)v  +  Ypp  +  (Yr  +  YrFF  +  YrFFF2  )r 
+  YV|V|  v  Ivl  +Yv|r,  v  Irl  +  Yrlr,  r  Irl  +Yr|v,  r  Ivl 
+  (  Y0  +  YiJ)F  F  +  YitiFF  F2  )  0 

+  (  YwC  v  +  Yv60  0  )  v  0  +  Yw4  v  r  0  +  (  Yrr0  r  +  Y^  0  )  r  0 
+  (  y8  +  Y888  52  +  Y8FF  F2  )  6  +  y5| r  1  6  I  r 1  +  ^  YFk8Fk 

+  (  Ycce  +  Yceeee  ^  +  YcceO0  •t*2  )  6 

Zcw  =  Zy  v  +  Zy,w  +  ZpP  +  Z^q  +  Z^.r  +  Zuu  +  Zvv  +  Zww  +  Zqq  +  ZI.r 

+  zz  zo  +  Ze  e  +  z5  8  +  ^  zFk8Fk 
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KCW=  (Ky+KyUu  +K^uuu2)v  +  Kwvv  +  KpjS  +  K^q  +  Ki.r 

+  (  Kv  +  KvF  F  +  KvFF  F2  )  v  +  Kp  p  +  (  Kr  +  KrF  F  +  KrFF  F2 )  r 
4-Kvivi  V  I v I  +  Ky | r |  v  Irl  +  Kr | r|  r  Irl  +Kr(v|  r  Ivl 
+  (  +  K^p  F  +  K^pp  F2  )  <{> 

+  (  KW0  v  +  KV00  <)> )  v  <}>  +  Kvr0  v  r  4>  +  (  Krr0  r  +  Kr66  <J> )  r  <(> 

+  (  K5  +  K555  82  +  K5Fp  F2 )  8  +  ^  Kpk8pk 

+  ( Kcce+  Kcceee  e2  +  4^  )  c2  e 

M<2\v=  My  v  +  M^w  +  Mpp  +  Mqq  +  Mj.r  +  Mu  u  +  Mvv  +  Mww  +Mqq  +  Mrr 
+  Mz  zq  +  Me  0  +  Mg  8  +  Mpk8pk 

NCW=  (Ny+N^u  +NVuuu2)v  +  NyyW  +  Npp  +  Nqq  +  lNi.  +  Nhlu  +  Nhillu2)r 
+  ( Nv  +  Nvp F  +  NvFp F2  )  v  +  Np p  +(Nr  +  NrpF  +  NrFFF2)r 
+  NV|V|  v  Ivl  +Nv|r|  v  Irl  +Nr|r|  r  Irl  +  Nr|v!  r  Ivl 
+  (  +  N^p  F  +  N^pp  F2 )  <J> 

+  (  Nw^  v  +  4> )  v  <)>  +  Ny^  v  r  <f>  +  (  r  +  <(> )  r  * 

+  ( Ng  +  Nggg  S2  +  Ngpp  F2 )  8  +  Ng| r  |  81  r  I  +  ^  Npk8pk 

+  ( Ncce+  Ncceee  e2  +  tp )  c2  e 

In  the  above  F  is  the  Froude  number,  defined  by  -s=  where  U  =  V  u2  +  v2  +  w2  is  the 

VgL 

instantaneous  forward  speed  of  the  ship  and  8  is  the  rudder  angle. 

Np  is  the  number  of  fins  installed,  8pk  is  the  angle  of  the  k1*1  fin.  The  fin  model  is 

included  to  account  for  the  effects  of  fins,  particularly  on  SWATHs.  As  currently  implemented, 
the  fin  angle  is  assumed  to  be  a  constant  which  is  set  at  the  beginning  of  the  run;  that  is,  no 
active  control  of  fins  is  considered.  The  fins  are  assumed  to  be  simple  lifting  surfaces  mounted  at 
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an  arbitrary  location  and  angle.  The  X  force  is  proportional  to  the  square  of  deflection,  5p;.  It  is 
assumed  that  Xp  is  minimum  at  8p  =  0. 

The  term  XRp  and  the  terms  with  coefficients  with  c  and  e  for  subscripts  can  be  used  to 
represent  the  resistance  and  propulsion  model.  Details  of  the  resistance  and  propulsion  model 
are  given  in  the  following  section. 

The  terms  included  in  this  hydrodynamic  model  were  selected  by  considering  the 
physics  of  the  flow  around  the  hull,  propellers,  and  rudders,  and  the  likely  functional 
dependence  of  the  resulting  flow  forces  on  the  mean  flow  over  the  rudder,  c,  and  the  effective 
angle  of  attack,  e.  Moderate  order  polynomials  in  these  variables  and  absolute  values  of  v  and 
r  were  considered.  Terms  were  included  in  the  model  only  if  the  resulting  contribution  would 
have  the  expected  symmetry  or  anti-symmetry.  That  is,  surge  force  is  symmetric;  sway  force, 
yaw  moment,  and  roll  moment  are  asymmetric.  Not  all  of  the  coefficients  in  the  model  are 
required  for  all  ships.  In  the  current  model,  some  asymmetry,  as  for  single-screw  ships,  may  be 
excluded.  In  this  case,  appropriate  terms  may  be  added  to  the  model. 

Some  terms  included  in  several  of  the  previous  models  mentioned  in  the  introduction 
would  not  satisfy  this  criterion.  Some  other  terms  are  ones  for  which  data  is  not  likely  to  be 
available.  The  appendix  lists  additional  terms  which  have  been  utilized  in  previous  models. 
These  terms  have  been  retained  in  the  computer  program  to  provide  flexibility. 

The  numerical  values  may  be  obtained  either  from  model  experiments  or  by  prediction. 

Smitt  and  Chislett[3],  and  reference  10  discuss  model  testing  techniques.  Acceleration  terms 
(Yy  ,  etc.)  may  be  predicted  with  reasonable  accuracy  by  potential  theory.  For  other 

coefficients,  predictions  are  usually  not  too  reliable  unless  a  similar  hull  is  available. 
Prediction  techniques  for  the  coefficients  are  beyond  the  scope  of  the  report  here.  For  example, 
see  reference  10  for  information  on  the  topic. 

Resistance  and  Propulsion  Model 

The  resistance  and  propulsion  model  calculates  the  net  surge  force  XRp  due  to  the  calm 
water  resistance  and  the  propellers.  In  this  section,  the  current  model  for  calculating  XRp  as 
well  as  two  older  methods  of  calculating  XRp  will  be  presented.  The  older  methods  are 

included  for  compatibility  with  older  sets  of  coefficients. 

The  current  method  of  calculating  XRP  has  several  advantages  over  these  older 

methods.  First,  flow  due  to  the  propellers  is  calculated  and  its  effect  on  the  rudder  forces  is 
modelled  explicitly,  making  possible  greater  accuracy  over  widely  varying  operating 
conditions.  Secondly,  effects  of  speed  on  the  non-dimensional  coefficients  should  be  eliminated, 
except  at  high  Froude  numbers,  and  thus  a  single  set  of  coefficients  should  suffice  for  all  speeds. 
In  contrast,  with  the  other  models  it  is  common  to  use  separate  sets  of  coefficients  at  different 
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speeds,  which  is  necessary  because  in  current  practice  the  resistance  and  propulsion  effects  are 
imbedded  in  coefficients  which  do  not  explicitly  depend  upon  Froude  number  (for  resistance) 
and  shaft  speed  (for  thrust). 

The  net  surge  force  due  to  calm  water  resistance  and  propulsion  is  given  by: 

Xrp  =  Np  Xqjrop  (1  - 1)  T  -  Ry 

where  Np  is  the  number  of  propellers,  t  is  the  thrust  deduction,  T  is  the  open  water  propeller 
thrust  per  shaft,  Ry  is  the  calm  water  resistance,  and  xCpTOp  is  a  coefficient  which  corrects  for 
possible  inaccuracies  in  the  calculated  values  of  T,  which  is  not  necessarily  well  known. 
Resistance  is  obtained  by  linear  interpolation  in  a  table  of  resistance  coefficients  which  are 
defined  by: 


Cd  =  ' 


Ry 


§U2L2 

Propeller  open  water  thrust  and  torque  coefficients  Ky  and  Kq  are  represented  in  the  form: 


T 

*T  =  =  ckt0  +  Jr  +ckt2  h2 


kQ  = 


(15) 


(16) 


p^5  =  <*0)+  <*Qj  Jq  +KQ2  Iq 
where  Jy  is  the  advance  coefficient  Ua  /  n  D,  n  is  the  shaft  speed,  and  D  is  the  propeller 
diameter.  The  speed  of  advance  Ua  is  given  by  U  (1  -  wt)  where  wt  is  the  thrust  wake  fraction. 
Jq  is  defined  similarly  for  torque  using  Wq,  the  torque  wake  fraction.  Given  U,  t,  w,  and  a 
machinery  model  as  discussed  in  the  following  section,  equations  (15)  and  (16)  can  be  solved  for 
nandT. 

Finally,  in  this  model  two  other  quantities  used  in  the  force  and  moment  model  of  the 
previous  section  are  also  calculated.  The  mean  flow  over  a  rudder  is  given  by: 


•■Vi 


Ap 

(1  -  wt)  u  +  k  UAtJ2  +  (  1  -  ^ )  (l'Wt)2  u2 


where  Ap  is  the  area  of  the  rudder  directly  behind  the  propeller  and  Ap  is  the  area  of  the 

rudder,  wt  is  the  wake  fraction  of  the  propeller  jet  velocity  at  infinity  and  k  =  UA  /Ua~> *s  a 
(xr-xp) 

function  of  where  x,.  is  the  rudder  location  and  *  3  the  propeller  location,  k  is  derived 

theoretically  in  reference  11.  Reference  7  contains  a  convenient  plot  of  this  function.  Ua<»  is 
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given  by: 


Uaoo  =  -O  -  wt) 


.  .  8Kr 

(1  -  wt)2u2  + - (n  D)2 


Finally,  the  effective  angle  of  attack  of  the  flow  over  the  rudder  is  given  by: 


e  =  5  -  tan"1 


v  r  xr 
c  c 


This  is  the  version  of  propeller-rudder  interaction  model  proposed  by  Abkowitz[7]  following 
Norrbin[6],  Use  of  c  and  e  makes  possible  a  more  accurate  calculation  of  the  influence  of 
propeller  speed  on  the  forces  and  moments  on  the  ship. 

The  first  of  the  two  older  models  is: 


Xrp  =  (1  - 1)  Np  T  -  Rt 


where  T  is  defined  by  equation  (15).  Rj  for  this  model  is  given  in  pounds  by: 


RT  = 


550  PE 
U 


where  PE  is  the  ship  effective  horsepower  represented  by  a  polynomial. 


The  propeller  thrust  is  calculated  as  in  the  first  model. 

The  second  of  the  two  older  models  utilizes  the  simple  expression: 

XRP  =  Xy  u  +  Xyj  T?  +  XM  "\p 


where  u  =  u  -  Uo  is  the  speed  loss  and  Ug  is  the  initial  forward  speed  of  the  ship.  In  this 
model,  hull  resistance,  propeller,  and  (implicitly)  machinery  dynamics  are  all  combined  into  a 
simple  function  of  the  speed  loss. 

Machineiy  Model 

The  machinery  model  governs  the  variation  of  propeller  speed  throughout  the 
maneuver.  There  exist  three  simple  models  and  two  detailed  models.  For  the  simple 
machinery  models,  either  speed,  torque,  or  power  is  held  constant.  The  detailed  machinery 
models  are  for  either  a  diesel  engine  or  a  turbine  and  include  the  effects  of  engine  dynamics. 
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governors,  torque  limiters,  and  drive  train  dynamics.  The  detailed  models  were  developed  by 
Propulsion  Dynamics,  Inc.[4]  and  will  not  be  included  here. 


WIND  FORCES 

The  wind  forces  and  moments,  Fy/IND/  are  ^  id  on  data,  Cy/uqp,  measured  in  a  wind 
tunnel  for  a  range  of  headings.  The  wind  force  and  moment  coefficients,  C^vind-  are  defined  by 

Pa  t  ->  Pa  i 

non-dimensionalizing  the  measured  forces  by  —  L2  U2  and  the  moments  by  —  L3  U2  where  pA  is 
the  mass  density  of  air. 

This  non-dimensional  data,  Cy/IND-  are  interpolated  linearly  in  heading  to  obtain  the 
values  for  the  correct  relative  wind  direction,  tyyy  -  \y,  where  yvv  is  the  wind  heading: 


VRW  =  tan'1 


Uw  sinfyvv  -  y)  -  vj 
Uw  cos(\yw-y)  -  uj 


(17) 


where  ui  and  vj  are  the  ship  velocities  in  the  inertial  coordinate  system.  FwiND  is  then 
calculated  by  dimensionalizing  the  resultant  values,  CwiND^VRw)-  The  force  components  are 


dimensionalized  by  ^  L2  Drw  and  the  moment  components  by  ^  L3  Ur^V  where 


Urw  =  V  l  UW  cos(y  w  -  V)  -  ui  ]2  +  [  Uw  sin (\yw  -  V)  -  vi  ]2  (18) 

The  wind  direction,  Vvv  is  measured  relative  to  the  xg-axis,  and  is  such  that  with  the  ship 
parallel  to  the  x^-axis,  Vvv  =  0°  is  a  following  wind,  Vw  =  90°  is  a  port  beam  wind,  and 
vyyy  =  180°  is  a  head  wind.  In  these  equations,  Uw  is  the  wind  speed  which  will  be  defined  in 
the  next  section.  These  calculations  are  performed  in  the  yawed  coordinate  system. 

Unsteady  Wind  Model 

An  unsteady  wind  velocity  model  is  incorporated  using  the  Davenport[12]  spectrum  for 
variation  of  the  longitudinal  component  of  the  wind.  No  lateral  spectrum  is  used  at  present,  so 
the  direction  of  wind  is  assumed  to  be  constant. 

The  unsteady  wind  velocity,  Uyy(t),  used  in  equations  (17)  and  (18)  is  given  by: 

Uw(t)  =  UwiND  +  V2  SWIND  <<“j)  AWj  cos(ci)jt  +  Jjj)  (19) 


where  UwiND  *s  the  mean  wind  velocity  in  knots  at  an  elevation  of  10  meters  above  the  water 
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surface  and  the  Davenport  wind  spectrum  density  function  Swind  ^  in  m2  /  sec  is  given  by  : 
1.459  x  105  k  to ; 

SwiND  toj)  =  (1_+  x^)4/3  ' 


(Oj  is  the  frequency  in  radians  per  second,  is  the  random  phase,  distributed  uniformly  over 
(0,2rt],  k  is  a  drag  coefficient  taken  as  0.005,  and 

x  =  1.378  xlO5 - — - 

UWIND 

In  the  model,  the  spacing  of  the  discrete  C0j  is  selected  so  that  each  component  in  the 
summation  of  equation(19)  has  the  same  amplitude;  the  frequencies  are  spaced  more  closely 
where  the  spectrum  is  larger. 


Approximation  for  Wind  Force  Coefficients 

If  measured  Cwind  ^ata  are  not  available,  an  alternative  is  the  model  applied 
previously  by  McCreight[5)  for  monohulls: 

XWIND  -  Y  Urw2  AF  Cxwmax  sin^|(  1  Vrw  1  *f  ^  1  Vrw  1  <n 


PA 


YWIND  =  Y  Urw2  aS  CYW  sin(  yRW  ) 

^  ¥Vmav 


max 


kWIND  =  Y  Urw2  As  Hs  CKw  sin(vRW)  cos2^ 


max 


Pa  l 

nWIND  =  —  URW  AF  CNw  sin(  2^rw  > 


max 


where  Ap,  As,  Hs,  Pa;  p,  and  L  are  projected  frontal  and  side  above  water  areas,  the  vertical 
coordinate  of  the  centroid  of  the  projected  side  area  of  the  ship,  air  mass  density,  water  mass 
density,  and  ship  length,  respectively. 

The  values  0.70,  0.80,  1.30,  0.10  have  been  adopted  for  the  wind  force  coefficients 

Cx  ,  Cy  /  Cr  ,  and  ,  respectively.  This  approximation  is  based  on 

wmax  wmax  wmax  wmax 

experimental  results  of  Aage[13]  for  a  number  of  ships.  The  roll  dependence  is  adapted  from 
Sarchin  and  Goldberg[14], 
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WAVE  FORCES 

A  seaway  acting  on  a  ship  causes  forces,  some  of  which  can  be  neglected  in  the 
maneuvering  problem  and  some  which  must  be  accounted  for.  The  incident  seas  can  be 
represented  using  the  model  of  St.  Denis  and  Pierson[15],  in  which  the  sea  surface  C(xo,yo-t) is 

represented  as  the  sum  of  sine  waves  of  varying  lengths  and  directions. 

C(x0,y0,t)  =  Re  jLj  e‘iKk(x0  cospj+yg  sinjXj )  +  i(Oj,.t  +  iejk  A|ij  At0k 


2 

“k 

for  short-crested  seas,  where  is  the  wave  number  —  of  the  kth  frequency  of  m 

S 

frequencies,  pj  is  the  jth  of  n  directions,  is  the  random  phase  of  the  wave  component  of 
frequency  coj.  and  direction  pj,  uniformly  distributed  over  (0,  2jc]  and  S^  is  the  Bretschneider[16] 
wave  spectrum  with  cosine-squared  spreading  given  by 

“  cos2  <4 '  Pc>  S^fo^)  0  <  I  p  -  Me 1  -  2 

=  0  Iclp-Mc1^" 

and 

-4  A 

_  ,  s  x2  A  c  1944.5Tnff-4 

S^<«»>k)  =  487.1  (Cw)  1/3  Tq  (o  5e  0 

2 

where  (£w)  3/3  is  the  significant  wave  height,  Tg  is  the  modal,  or  peak  energy,  period  and  pc  is 

the  (predominant)  wave  direction.  The  wave  direction  p  is  measured  relative  to  the  xg  axis, 
and  is  such  that  with  the  ship  parallel  to  the  xg-axis,  p  =  0°  is  a  following  wave,  p  =  90°  is  a 
port  beam  wave,  and  p  =  180°  is  a  head  wave.  This  convention  is  the  same  as  for  the  wind 
heading. 

For  long-crested  waves  the  appropriate  expression  is: 

C(x,y,t)  .  Re  £,  eiKK(x0  «>  W0(t)  SinUj>  ♦  tokt  *  iek  V2S?((])k)  4l0|< 


where  p  is  the  wave  direction  and  is  the  random  phase  uniformly  distributed  over  (0,2ji]. 

Ogilvie[17]  presents  an  excellent  survey  of  the  resulting  forces  and  methods  of 
predicting  them.  Here  we  just  give  a  brief  summary  of  these  forces  and  discuss  the  assumptions 
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and  approximations  used  in  the  program  for  calculating  the  slowly-varying  wave  drift  force 
components,  as  well  as  indicate  why  the  other  terms  were  omitted. 

For  purposes  of  explanation,  we  now  consider  two  plane  waves  of  different  lengths 
propagating  in  the  same  direction 

£(t)  =  A}  coso>|t  +A2  coso^t 

Assuming  a  perturbation  expansion  in  the  wave  amplitudes  Aj  and  A2  and  ignoring  terms 
greater  than  quadratic,  the  forces  are  of  the  form: 

F(t)  =  Re{A}  G](cO|)  e  ^+A2Gj((o2)  2  }  +  Re(A^  G2(ci)-],-(Oj)  +  A2  G2((o2,-cd2)  ) 

1  _  , .  2  _  ,  ,  i2(0i t  1  ,,2  _  ,  .  i2a>it , 

+  -  Re{A-,  G2(co1,co1)  e  )  +  2  ReiA2  G2(£o2,o)2)  e  L  } 

+  Re(Aj  A2G2(o>1,co2)  el(&>1  +  W2)t}  +  Re{A^  A2  G^oip-o^le  i(0>1  (20) 

In  equation  (20)  the  first  and  second  terms  which  are  proportional  to  A]  and  A2,  respectively, 
are  the  linear  exciting  forces,  which  will  be  ignored  in  the  present  work.  These  are  responsible 
for  the  motions  of  ships  in  waves  as  predicted  by  the  standard  seakeeping  programs,  but  the 
motions  are  periodic,  average  out  to  zero,  and  thus  are  not  of  interest  for  the  pure  maneuvering 
problem. 

The  remaining  terms  are  the  nonlinear  wave  forces.  These  are  quite  small  in  magnitude 

relative  to  the  first  order  forces,  and  are  generally  neglected  in  standard  seakeeping  studies. 

Several  of  these  are  important  for  stationkeeping  and  trackkeeping  studies,  because  while  they 

are  small,  their  mean  is  non-zero  and  can  cause  considerable  excursions  in  surge,  sway,  and  yaw, 

2 

where  there  are  no  restoring  forces  to  balance  them.  The  next  two  terms,  proportional  to  A-j  and 
2 

A2,  respectively,  are  the  mean  drift  forces,  and  are  important.  The  next  three  terms, 
proportional  to  e'20*1*,  e*20*2*,  and  e'^  +  °>2^t,  respectively,  act  at  even  higher  frequencies  than 

i(<Di  “  G>2)t 

the  incident  waves,  and  thus  are  negligible.  The  last  term,  proportional  to  A]  A2  e  , 

results  in  a  slowly- varying  drift  force  which  is  often  also  neglected,  but  may  be  important  in 
trackkeeping  with  strict  requirements  on  the  amount  of  drift  permitted  from  the  nominal 
trackline.  This  is  also  important  for  stationkeeping,  but  is  traditionally  ne6.txicu  in  the  static 
force  balance  approach  often  used  in  these  problems.  In  the  two  wave  case  considered  here,  the 
slowly-varying  force  is  a  simple  sinusoidal  beat  at  the  difference  frequency.  The  two  main 
expressions  for  forces  of  equation  (20)  can  be  extended  to  represent  the  responses  to  an  irregular 
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seaway  in  terms  of  single  and  double  integrals  (summations  for  the  discrete  case).  Most  of  the 
terms  required  in  such  an  expression  are  unknown,  however,  and  an  approximate  method  must  be 
found,  utilizing  only  the  mean  drift  force  operator  G2(w,-co). 

Based  on  the  discussion  given  above,  only  the  slowly-varying  drift  forces  will  be 
included  in  this  model.  That  is, 

Fwave  =  Fsvdf 

where  FsvDF  will  be  defined  in  the  next  section  utilizing  the  mean  drift  force  operator 
Co(a>,-co)  which  may  be  obtained  from  experimental  data  or  from  Dredictions. 

Slowly-Varying  Wave  Drift  Force  Model 

In  this  program,  slowly-varying  wave  drift  forces  in  short-crested  waves  are  calculated 
by  an  approximate  method  due  to  Marthinsen  [18].  This  method  is  based  on  the  assumption  of  a 
narrow-banded  spectrum  and  the  idea  of  a  local  wave  frequency  and  direction  to  enable  the  use 
of  regular  wave  drift  force  operators  as  the  basis  of  the  calculation.  This  is  a  more  systematic 
and  general  approach  which  is  equivalent  to  the  approximation  of  Hsu  and  Blenkem[19].  In 
the  case  of  a  long-crested  wave,  Marthinsen  showed  through  several  examples  that  this 
approach  is  equivalent  to  the  commonly  used  approximate  methods  of  Newman[20].  This 
approximation  uses  only  the  regular  wave  drift  force,  which  is  the  most  readily  available  from 
experiment  or  calculation.  These  forces  are  calculated  in  the  yawed  system. 

We  will  not  derive  the  method  here.  The  calculation  method  is  as  follows: 

Calculate  the  Hilbert  transform  of  the  wave  elevation,  which  is  given  by 

„(*0.y0,O  -  In.  j£,  g.,  ^  “>k  VS{tok.»j)41ij4®k 

for  the  short-crested  case  and  by 

„(x0,y0,t>  -  ta  g.,  e*‘‘(x°  COS4i*yo  Sin“i)  +  ‘“k' +  *'k  V2S;(ck)  40)k 


for  the  long-crested  case.  Then  the  instantaneous  wave  amplitude,  frequency,  wave  number 
components,  and  wave  direction  are  given  by: 


A(t)  =  V  C( t)2  +  Tj(t)2 


<o(t)  = 


M*>  = 


q(t)  gt 


n(t) 


sect) 

8x 


c«> 


8t 

Sq(t)  N 

8x  J 


/  A^t) 
/  A2(t) 
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The  resulting  approximate  expressions  for  the  slowly-varying  drift  forces  and  moment  are 
XsVDF  =  \  P  g  L  A2(t)  Cx(  oft).  y(t)  ) 

YsVDF  =  \  P  g  L  A2(t)  CY(  co(t),  v(/(t) ) 

NsVDF  =  \  P  g  L2  A2(t)  Cn(  (0(t),  vy( t)  ) 

where  in  the  previous  notation 

Go  ( to,  —(0  ) 

Cx(  co(t),  y(t)  )  = — - - 

2  pgL 

represents  the  non-dimensional  mean  surge  drift  force  for  a  regular  wave  of  frequency  w  and 
direction  vy  and  Cy  and  Cpj  are  similarly  the  non-dimensional  drift  force  and  moment  for  sway 

and  yaw,  respectively. 

In  the  simulation,  the  drift  forces  are  updated  only  on^e  per  second  and  are  assumed 
constant  during  the  one  second  time  interval.  Further,  in  the  initial  phase  of  the  simulation, 
the  forces  are  ramped  up  from  0  at  t  =  0  to  the  full  value  at  t  =  t^^p,  which  is  an  input 

variable. 

Strictly  speaking,  the  slowly-varying  drift  force  model  is  valid  only  for  zero  forward 
speed,  because  experimental  data  is  available  only  at  zero  speed.  Utilizing  the  model  at  low 
speeds  is  reasonable,  however.  At  high  speeds,  drift  forces  are  less  significant  because  the  hull 
forces  which  are  proportional  to  speed  squared  are  dominant  and  can  easily  compensate  for  the 
effect  of  drift  forces. 

TOW  FORCES 

The  force,  Fj0w>  due  t0  towing  a  body  through  the  water  is  also  modelled.  It  is 

assumed  that  the  tow  cable  is  parallel  to  the  total  ship  velocity,  U,  and  that  the  force  is 
proportional  to  the  square  of  the  total  ship  velocity.  The  tow  point  is  located  at  (xj,  yj,  zj). 

kTOW  =  TyzT 


MTOW  =  Tx  ZT 


NTOW  =  TyxI  -  TxYT 


where 


Tx  =  cos  [3  Djow 


2 

UjOW 


Tv  =  sin  P  Djqw 


U2 


>'  2 

U70W 


and  Djow  1S  the  drag  at  a  reference  speed  of  UjqW' 

SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  program  first  initializes  the  state  vectors  for  the  ship  and  then  the  equations  of 
motion  are  integrated  forward  in  time  to  obtain  the  specific  solution  for  a  given  set  of  initial 
conditions  and  rudder  and  speed  commands.  Turns,  spirals,  zig-zags,  and  arbitrary  sequences  of 
rudder  and  speed  changes  may  be  commanded. 

The  initial  conditions  for  the  time  domain  solution  are  chosen  to  represent  the  steady 
state  solution  corresponding  to  the  ship  moving  at  constant  speed  and  heading,  subject  to  the 
wind  and  waves  assumed.  Mean  values  of  unsteady  wind  and  slowly-varying  wave  drift  forces 
are  used  during  the  initialization.  The  unsteady  parts  are  ramped  up  from  zero  to  full  value  at 
a  (user-specified)  time  after  the  beginning  of  the  run.  This  is  necessary  in  order  to  prevent 
spurious  starting  transients  in  the  solution 

The  initial  conditions  are  determined  by  minimizing  the  sum  of  the  squares  of  the 
accelerations  in  six  degrees  of  freedom.  This  yields  rudder  angle,  shaft  horsepower,  heave, 
roll,  pitch,  and  yaw  (drift  angle)  to  balance  wind  forces  and  forces  due  to  asymmetry  of  the  ship 
calm  water  hydrodynamic  model. 

The  method  adopted  for  finding  initial  conditions  yields  solutions  which  result  in 
negligible  starting  transients.  This  has  been  demonstrated  by  running  calm  water  solutions  with 
no  steering  or  speed  commands  for  periods  of  time  up  to  1000  seconds  (full  scale)  which  resulted 
in  essentially  no  change  in  the  ship  speed,  heading,  or  attitude. 

At  any  time  t,  the  equations  of  motion,  see  equations  (1-6),  can  be  solved  algebraically 
for  the  acceleration  in  terms  of  the  velocities,  displacements,  and  forcing  functions.  Force  and 
moment  components  for  hydrostatic,  wind,  wave  drift,  and  tow  forces  are  transformed  to  the 
maneuvering  coordinate  system  using  equation  (13)  and  then  added  to  the  calm  water 
hydrodynamic  forces  as  in  equation  (14).  Together  with  the  kinematic  relations  of 
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equations  (7-12)  for  the  time  derivatives  of  the  position  and  orientation,  and  for  certain 
variables  used  in  the  (optional)  machinery  models,  these  form  a  system  of  differential 
equations  which  are  integrated  numerically  from  the  initial  conditions  as  described  in  the 
previous  section  to  find  the  motion  of  the  ship  as  a  function  of  time.  In  the  present  program,  an 
adaptive  Runge-Kutta  method  developed  by  Merson[21]  is  used. 

Test  runs  in  calm  water  conditions  showed  that  predictions  of  sharp  turns  (rudder  angle 
35  degrees)  were  practically  independent  of  step  siz.  for  step  sizes  ranging  from  0.0333  seconds 
to  5.0  seconds.  For  step  sizes  greater  than  1.0  seconds,  the  step  reduction  feature  came  into  play 
in  the  initial  portion  of  the  run  when  accelerations  were  largest. 

CONCLUSION 

This  report  describes  a  mathematical  model  to  predict  the  motion  of  a  surface  ship 
while  maneuvering  in  calm  water,  including  unsteady  wind  and  slowly-varying  wave  drift 
forces,  and  forces  due  to  a  towed  body.  This  model  incorporates  a  model  of  propeller-rudder 
interaction  based  on  the  physics  of  the  situation,  rather  than  a  simple  curve  fit  to  experimental 
data  as  has  been  used  previously.  This  mathematical  model  has  been  implemented  in  a 
computer  program,  MPSS  (Maneuvering  Program  for  Surface  Ships).  Hydrodynamic  coefficients 
must  be  provided  for  the  mathematical  model.  In  addition,  terms  used  in  previous  programs 
which  are  not  included  in  the  model  here  may  optionally  be  included.  This  program  is 
applicable  to  SWATH  ships  as  well  as  monohulls. 
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APPENDIX  -  ADDITIONAL  TERMS  IN  WAVMAN44  OR  TOWMAN87 

The  purpose  of  this  appendix  is  to  present  terms  which  are  included  in  the  calm  water 
forces  and  moments  in  WAVMAN44  or  TOWMAN87,  but  are  not  included  in  the  model 
presented  in  this  report.  An  option  exists  in  MPSS  so  that  these  additional  terms  can  be 
utilized.  This  makes  it  possible  to  run  MPSS  using  old  data  sets. 

VVAVMAN44  is  a  four-degrees-of-freedom  program  so  that  heave  and  pitch  are  not 
included.  In  WAVMAN44,  U  is  defined  as  V  u2  +  (w  cos  b)2  ,  rather  thanV  u2  +  v2  +  w2.  This 
U  is  used  to  dimensionalize  the  forces  and  moments.  In  addition,  WAVMAN44  includes  some  of 
the  terms  which  are  multiplied  by  F  in  the  model  presented  in  this  report.  However  in 
WAVMAN44  those  terms  are  multiplied  by  u.  ii  is  defined  as  (u  -  Uo  )  and  Ug  is  the  initial 
forward  speed  of  the  ship.  When  the  WAVMAN44  option  is  utilized  in  MPSS,  the  definition 

of  U  stated  here  and  the  substitution  of  u  for  F  will  be  implemented. 

WM  WM  WM  WM 

The  additional  terms  in  WAVMAN44,  given  as  ,  Ycw,  K(~.y,  and  Ncw  are: 

WM  ,  t 

XCW  =Xvv  +Xrlvl  r  1  v I  +XVVFPV2F2 

>Xu55F^  52p  +  x8  8  +  X55u  &  u  +  X6uu6°2 

+  Xv5v8  +  Xr6r5  +  X|5|v  181  v  +  Xg  |  v  |  8  Ivl  ♦  Xg^v2 

YCW  =Y0  +  Y0uii  +  Yvvv2  +  Yvvvv3  +  YVITvr2  +  Yrrrr3 

+  Yrvvrv2  +  Yry/Tvi  rVTv’l  +  )r|  4>  Iri  +  Y|0|r  l<})l  r  +  Y^  b3  +  Ygp  Sb 

+  Y^g  v  r  8  +  Y5vv  8  v2  +  Y6rr  5  r2  +  Ygg  S2  +  Yvgg  v  S2  +  Yrgg  r  S2  +  Y5vt  8  v  r 

+  Y6|V|  8  Ivl  +YV|5|VI  8|  +  Y5a8u  + Y55a82u  + Y656uS3u 

^■CW  =  ^pu  P  u  +  Xp|  pi  p  Ipl  +Kp|piuP  Ipl  ^  +  ^rrr  ^  +  ^rrv  r  v 

+  K ^  r  v2  +  KvF  v  F  +  KvFF  vF2  +  Kp/p  p  /  F  +  KpF  p  F  +  Kppp  p  F2 

+  K,(>|r  Ibl  r  +  K^ij-i  b  Irl  +  b’  +  Kgiy,  8  Ivl  +  Kg^u 
WM  a 

^CW  =  No  +  Nqu  +  Nw  v  +  Nyw  v3  +  Nvrr  v  ' 
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+  Njyj.  r2  +  Nj^r  v^N^VhT.  +N^(r|  0  I r I  +  N|*|r  f 0 1  r+N^^ 

+  N0i5|  0  181  +  Nvr5vr8  +  N5vv8v2  +  N5n.8r2  +  Nss82 

+  Nv5Sv82  +  Nr55r82  +  N5vr8 vr  +  NSiv|  8  Ivl  +Nv,6|vl  8|  +  N6G8u 

+  N5SuS2u  +  N5S5g83u 

While  TOWMAN87  includes  some  heave  and  pitch  calm  water  hydrodynamic  and 
hydrostatic  terms,  it  does  not  include  the  necessary  inertial  terms  and  transformations  between 
coordinate  systems  required  for  a  six  degrees  of  freedom  model. 

TOWMAN87  includes  calm  water  forces  for  heave  and  pitch.  However,  it  does  not 
include  the  inertial  terms  and  coordinate  transformations  required  for  a  consistent  six  degrees  of 
freedom  model.  In  TOWMAN87,  U  is  defined  as  Vu2  +  (w  cos  0)2  ,  rather  thanV u2  +  v2  +  w2 . 
This  U  is  used  to  dimensionalize  the  forces  and  moments.  In  addition,  TOWMAN87  includes 
the  terms  which  are  multiplied  by  F  in  the  model  presented  in  this  report.  However  in 
TOWMAN87  those  terms  are  multiplied  by  powers  of  U.  When  the  TOWMAN87  option  is 
utilized  in  MPSS,  the  definition  of  U  stated  here  and  the  substitution  of  U  for  F  will  be 
implemented  when  the  flag  indicating  TOWMAN87  data  is  set  in  the  input.  The  additional 
terms  in  TOWMAN87,  given  as  ^  Y^,  Mcw,  and  Ncvv  are: 


^CW  -  ^8  S  +  X8U  8  U 

YCW  =  Yvw  v3  +  Yv0  v  Ypj,  r  0  +  Y5v  8  v 
TM 

KCW  =  ^v0  v  ^  +^T0  r<^ 

TM 

MCW  =  <t> 

TTvI 

Ncw  =  Nrv  r  v  +  Nvw  v3  +  Nv0  v  Nr0  V  0+  N6v  8  v 


I 


% 
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